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Abstract 

We apply Frobenius integrability theorem in the search of invariants for one-dimen¬ 
sional Hamiltonian systems with a time-dependent potential. We obtain several 
classes of potential functions for which Frobenius theorem assures the existence of 
a two-dimensional foliation to which the motion is constrained. In particular, we 
derive a new inhnite class of potentials for which the motion is assurately restricted 
to a two-dimensional foliation. In some cases, Frobenius theorem allows the explicit 
construction of an associated invariant. It is proven the inverse result that, if an 
invariant is known, then it always can be furnished by Frobenius theorem. 
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1 Introduction 


The search for invariants (also called constants of motion or first integrals) for 
dynamical systems is a classical topic in nonlinear science, and in the past two 
decades there has been an intensive research on the subject. A comprehensive 
review of the recent advances on the held can be hnd in [^. There are several 
techniques for the construction of invariants, among which the direct method |Q, 
application of Noether’s theorem Q, utilization of Lie symmetries [^, Painleve’s 
analysis P] and Darboux method 0. None of these methods has an universal 
character, and in most cases one or more ad hoc assumption has to be made for 
obtaining concrete results. For instance, in the case of Noether and Lie’s approach 
an useful simplihcation is achieved by considering point symmetries only 0, |^. 


In the present work, we focus our attention on Hamiltonian systems with a 
time-dependent Hamiltonian function of the form 

= ^p'^+ V{q,t), (1) 


where V (g, t) is a time-dependent potential function and p and q are coordinates 
on a two-dimensional phase space. Hamiltonians of type (® appear in several 
branches of physics such as plasma physics and quantum mechanics. Here, we 
are interested in particular in the construction of invariants for suitable classes of 
potentials. By invariant, we understand precisely a function / = I{q,p,t) globally 
dehned and having the property 


dl dl dV dl 
dt ^ dq dq dp 


( 2 ) 


where the dot represents time differentiation. In other words, equation (|^) says 
that the function I is constant along the trajectories of the canonical equations 
of motion. The existence of an invariant immediately prevents the appearance of 
chaos, since when a hrst integral is available the motion is restricted to a two- 
dimensional surface I{q,p,t) = cte. As is well known, if the independent variable 
is continuous, chaos may take place only in dynamical systems with three or more 
dimensions. Note that an explicit time-dependence of I does not modify at all the 
fact that chaos is impossible if an invariant is known to exist (provided I{q, p, t) = 
cte. dehnes a sufficiently smoth surface). 


To begin our analysis, we write the dynamical vector held associated to the 
canonical equations of motion. 


d d dV d 
dt ^ dq dq dp 


(3) 
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Now, the definition reads u(/) = 0. The dynamical vector field u is defined 
in a three-dimensional space, with coordinates g, p and t. Incidentally, some years 
it was introduced a method for the derivation of invariants for three- 
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ago 

dimensional dynamical systems. This method relies on Frobenius integrability the¬ 
orem, and has been used in the analysis of the three-dimensional Lotka-Volterra 

0 ] 


system and the May-Leonard and Lorenz systems 


Presently, we apply 


Frobenius method in order to derive integrable classes of Hamiltonians of the form 
i). The large number of integrable three-dimensional Lotka-Volterra systems de¬ 
tected by Frobenius method is a good reason for optimism. Indeed, as seen in 
the following, we are able to find several integrable potentials by using Frobe¬ 
nius method with some auxiliary hypothesis. In particular, we find a new whole 
infinite family of potentials for which the motion is shown to be restricted to a 
two-dimensional surface. Finally, we prove the result that, if an invariant is known 
to exist, then it can be derived by Frobenius method. This last result shows that 
the Frobenius approach for invariants for one-dimensional time-dependent Hamil¬ 
tonian systems necessarily reproduces all known integrable cases in the literature. 

The article is organized as follows. In Section H, we describe Frobenius method 
and derive the basic equation to be solved in the case of Hamiltonian systems with 
Hamiltonian of the form (||). In Section HI, we use a linear ansatz in momentum 
for obtaining particular solutions of the basic equation derived in Section H. As a 
result, we derive the invariant linear in momentum for the forced time-dependent 
harmonic oscillator and Sarlet’s invariant |^-|T^. In Section IV, we con¬ 
sider functions rational in momentum, so obtaining an infinite family of potentials 
amenable to Frobenius method. As particular cases in this family, we found the 
potentials having an invariant quadratic in momentum and Giacomini poten¬ 
tials |]T^ |T^. Also we show a class of “weakly integrable” potentials, which have 
not in general an invariant but for which Frobenius method works. Section V is 
reserved to the conclusions. 


2 Frobenius method and basic equation 


The Frobenius method was already exposed in [^. Here, we only show the es¬ 
sentials of the procedure, with a view to application to one-dimensional time- 
dependent Hamiltonian systems. The basic aim in the Frobenius procedure for the 
construction of invariants of motion is to find a vector field v given by 


V = 


,9 

+ B — 
dt dq 


C 


dp 


(4) 


which is compatible with the dynamical vector field. In equation (^, A = A{q, p, t), 
B = B{q, p, t) and C = C{q, p, t) are functions to be determined. Compatibility 
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means that 


[u, v] = au + /3v , (5) 

where [, ] represents the Lie bracket, [u, v] = uv — vu, and a = a{q,p,t) and 
(3 = j3{q,p,t) are functions (not necessarily constants) on phase space and time. 
Equation (§) implies that the dynamical vector held, the compatible vector held 
and their Lie bracket are linearly dependent. If /? = 0, the compatibility condi¬ 
tion reduces to the statement that v is a generator of Lie symmetries for the 
dynamical system |^. For /? 7 ^ 0, there is a generalization. 

If a compatible vector held can be found, then the Frobenius theorem assures 
the existence of a maximal connected leaf which is tangent to both u and v, in 
each point in the three-dimensional space with coordinates (g, p, t). Moreover, the 
union of all such leaves provides a foliation, that is, each point {q,p, t) is contained 
in one and only one leaf. However, this two-dimensional foliation is not necessarily 
dehned in terms of the level surfaces of some function I{q,p,t). Nevertheless, for 
the derivation of an invariant common to the dynamical and the compatible vector 
helds we assume the existence of such function, so that I{q,p,t) = cte. provides a 
representation of the foliation. The fact that u and v are tangent to the maximal 
connected leaf implies that u(/) = v(J) = 0 . 

Suppose that a compatible vector held v is available. The natural question 
in this circumstance is how to construct an associated invariant I{q, p, t). The 
procedure for the derivation of an invariant in Frobenius method is identical to 
the construction of hrst integrals by Lie symmetry methods |@. By dehnition, 
u(/) = v(/) = 0 , since the maximal surface is tangent to both the dynamical and 
the compatible vector held. Consider 

v(/) = 0. (6) 

Denoting the integral surfaces of the compatible vector held by f{q,p,t) = cte., 
g{q,p,t) = cte., so that v(/) = v(p) = 0, we have that 

I = Hf,9) (7) 

is the general solution to (H). As shown in 0, substitution of the form in the 
condition u(/) = 0 leads to a hrst-order ordinary diherential equation, 

| = A(/.s). (8) 

where A is a function depending on / and g only. Expressing the constant of 
integration of the later equation in terms of f{q,p,t) and g{q,p,t), we obtain the 
invariant associated to the compatible vector held. As seen by the reasoning in 
this paragraph, the Lie and Frobenius strategies for derivation of invariants are 
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indeed equal, once v is known. The sole difference between the two methods is 
that the compatible vector field is not necessarily a generator of Lie symmetries. 

Equation is a first-order linear partial differential equation which has to 
be solved by the method of characteristics. One may expects, a priori, that this 
equation is more easy to handle than the original dynamical system. Otherwise, 
existence of a compatible vector field only assures that the motion is restricted to a 
two-dimensional foliation, but does not provides the explicit form of an associated 
invariant. Another possible drawback in the procedure is the resolution of (H). 
Systems having a compatible vector field but for which one is not able to find 
explicitly the associated constant of motion were termed weakly integrable in the 
literature [^]. In Section IV, we show an infinity family of potential functions which 
pertains to such weakly integrable category. 


Turning our attention to our specific problem, the study of one-dimensional 
time-dependent Hamiltonian systems, let us consider the dynamical and compati¬ 
ble vector fields (|^-^). Actually, there is no loss of generality in taking A = 0 and 
H = 1 in the definition of the compatible vector field. Indeed, suppose that the 
compatibility condition (j^) is fulfilled. Defining a new vector field v' by 


v' = /iU -|- z/v , 


we get, from (||), 

[u. v'] = (u(^) + 

Thus, u and v' are compatible. Now, define 

_ A _ 1 

^ Ap — H ’ B — Ap 


(9) 


( 10 ) 


( 11 ) 


We always take B ^ Ap, so that there is no singularity in the denominator of p 
and p as given in (^TJ). The case B = Ap can be shown to give only trivial results 
(dynamical and compatible vector fields linearly dependent). With the choice (|Tl|), 
we have 


f C + AdV/dq X ^ 
y B — Ap j dp 


( 12 ) 


For all compatible vector fields v, we can construct an associated compatible vector 
field v' as defined by (^). In conclusion, vector fields of the form 


V = 


A 

dq 


+ C{q,p,t) 


dp 


(13) 


are all that is needed in our calculations. There is no real gain in considering 
compatible vector fields more general than (1^. Observe that the use of com¬ 


patible vector fields like (p!3D is analogous to the use of generators in evolutionary 
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form, in the theory of Lie extended groups [Q. However, here there is even more 
simplihcation, since the coefficient oi d/d qm. \ was set to unity. Finally, a major 
advantage of using the reduced form (|^ is that the integral surface g = t = cte. 
of the compatible vector field is known in advance. There remains only the task 
of computing the second integral surface. 


Let us study the consequences of the compatibility condition. Computing the 


Lie bracket of the dynamical vector field and v as given in (p!3D , we find 

[u, v] = -C— + u(C') + —^ — . 


dq 


dq'^ j dp 


(14) 


Now using the compatibility condition (|^), we determine the coefficients a and (3, 

a = 0, (3=-C. (15) 


In general (3 is not zero, so that the compatible vector field is not a generator of 
Lie symmetries. 


The compatibility requirement (^) comprises three equations, since we are deal¬ 
ing with vector helds in three dimensions. Two of these equations were used to 
calculate a and (3. The third equation gives 


u(C') + 


d^V 

dq^ 


(16) 


The latter equation may be used to determine simultaneously the potential and 
the function C, and is the basic tool in Frobenius method for one-dimensional 
Hamiltonian systems with a time-dependent potential. In the next sections, we 
show several solutions for ([T^j. To deal with (|T^) we have to suppose appropri¬ 
ate particular forms for the function C. As remarked in the introduction, most 
methods of derivation of invariants are not free from some assumption. Frobenius 
approach is not an exception. 


To end this section, we comment on the possible equivalence between the meth¬ 
ods of Lie and Frobenius for the construction of invariants. Indeed, if a compatible 
vector field v is known, one may wonder if a generator V of Lie symmetries may 
be obtained. So, consider v' of the form (j^). Frobenius method generalizes Lie’s 
approach by allowing a non zero (3 function in equation (^. A look at equation 
0 ) shows that we have a Lie symmetry if we choose v so that 

u(z/) += 0 . (17) 

Any vector field v' of the form (^, with arbitrary /i and with p satisfying (O) , is a 
generator of Lie symmetries provided v is a compatible vector field. However, the 
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difficult point in this reasoning is that (^) is a first-order linear partial differential 
for u, which can be solved in general only if the integral surfaces of the dynamical 
vector field are known. Obviously, a priori this information is not available, at 
least for nontrivial systems. 

In Sections III and IV, we solve the fundamental equation (|T^) for selected 
functional forms of C{q,p, t). We use C in the form of functions linear and rational 
in momentum. These two cases are treated separately. 


3 C{q,p,t) linear in momentum 


As remarked earlier, to obtain definite results some assumption must be made on 
the functional dependence of C{q,p,t). In this section, we consider 

C = Co + C,p, (18) 


where Cq = Co{q, t) and C\ = Ci(g, t) are functions depending only on position and 
time. Inserting the linear ansatz ([T^) into ([T^) , there results that a quadratic poly¬ 
nomial in momentum must be identically zero. Equating to zero the coefficients 
of the various powers of p, we find a system of partial differential equations. 


dt 


+ 


dq 


dCo 

dt 


dq 

= 0, 

(19) 

+2C'oC'i 

7 

= 0, 

(20) 

dV d‘^V 
dq dq'^ 

= 0. 

(21) 


This system is nonlinear. To solve it, first observe that 
solution, 

6*1 = 0 , 


and 


6*1 = 


q-a 


admits two types of 

( 22 ) 

(23) 


where a = a{t) is an arbitrary function. The two branches will be studied 

separately. 
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3.1 Linear invariants 


Let us first consider the case Ci 
solved, yielding 


0. In this circumstance, (|^) may be readily 


6*0 = 


P 

P 


(24) 


where p = p{t) 7 ^ 0 is an arbitrary function. The form (p4D was chosen for later 
convenience. We can now solve for the potential using obtaining 


^ = Vo{t) 


F{t)q 

P 



(25) 


where Vo(^) and F{t) are new arbitrary functions of time. In fact, we can set 
Vo = 0 without loss of generality, since it does not contribute to the equations of 
motion. 


The potential (^5] ) corresponds to the forced time-dependent harmonic oscillator 
system. To construct an invariant for it using Frobenius method, consider the 
associated compatible vector held, which, according to ( 0 ) and (P!BD, reads 


d p d 
dq p dp 

The integral surfaces of v are specihed by 


(26) 


f =P -, 9 = t. 

P 


(27) 


Computing df /dt along the trajectories of the dynamical system, we obtain 


^ = -pl + L 

dg P P ' 


(28) 


As p = p{g), the right hand side of the last equation is indeed a function of / and 
g only, in accordance with the general result (§). Solving (|^) , we obtain the well 
known |]^ invariant linear in momentum for the forced time-dependent harmonic 
oscillator, 

I = pp — pq — F . (29) 


3.2 Sarlet’s potential 

The previous calculations show how the Frobenius method works. Let us consider 
the less trivial solution (^) for (|^ and proceed to the solution of (pnj-pil). Taking 
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into account equations 


and 
C^o = 


, we obtain 
a{q - cr) + 27 


{q - af ’ 

where 7 = 7 (t) is an arbitrary function of time. With Co and Ci specihed by 
and (|30D , we can hnd the potential using (| 2 l|). It reads 


(30) 


V = Vo(t) - a(q - cr) 


PQ 


7^ 


- 7 log(g - cr) , 


2p 2 (g-cr )2 

introducing the new arbitrary functions Vb = Vo(^) and p = pit) ^ 0 . 


(31) 


Having solved the system ( p!9| - |21|) , there remains the task of obtaining the asso¬ 
ciated invariant using the compatible vector held, which is in the present case 

^ (g-a)(p-d) -27 

^ dq {q — cr )2 dp 

The associated invariant surfaces are specihed by 

p — & 7 


(32) 


/ = 


g = t. 


q-a (g - cr )2 ’ 

Along the canonical equations of motion with potential (^), we obtain 


(33) 


df _ : P 

dg p ’ 

a Riccati equation. Note that the general result 
p{g). The general solution for (^) is 

1 




p ppi-j‘dt'/fP)’ 


(34) 

is indeed verihed, since p = 

(35) 


where I is the constant of integration for the Riccati equation. Actually, I is the 
invariant of the problem, and, using (|^) and (^5|) , we have 

j ^ f^_ (9 - ^)/P _ 

J p2 p(^p -a)- p{q -a)- 7p/ (q-a)' 


(36) 


The potential (^) and the invariant (pq) are not new, being derived by the hrst 


time by Sarlet and then by a variety of methods [12 


In this section, we have shown how Frobenius procedure works, in the elementary 
case of functions C{q,p,t) that are linear in momentum. In this way, we have 
obtained some already known results. In order to hnd new integrable or weakly 
integrable potentials, a more complicated momentum dependence of C (g, p, t) must 
be utilized. It can be easily shown that higher order polynomial forms of C{q,p, t) 
does not gives noting new. Hence, we proceed diherently in the next section, taking 
a rational form for C{q,p,t). 
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4 C{q,p,t) rational in momentum 


In this section we take the following ansatz for the solution of the basic equation 


( 0 ), 


_ Co + Cip 
P-C2 


(37) 


where Cq, Ci and C 2 are functions of coordinate and time only. We also assume 
that Co + C 1 C 2 7 ^ 0, so that dC/dp 7 ^ 0. Substitution of ( |37[) into (|^) yields an 
equation implying that a cubic polynomial in momentum is identically zero. Again, 
the coefficient of equal powers in momentum must be zero, resulting in a system 
of partial differential equations for Co, Ci, C 2 and the potential. Considering the 
coefficient of we get that Ci is a function of time only. For latter convenience, 
we represents the solution as 

Cl = p/p, (38) 


where p = p(t) is an arbitrary nonzero function of t. Taking into account (^), one 
can show that the term proportional to p^ in the basic system of equations yields 


Co =-d - ^(g - cr) - ^C 2 - ^ , (39) 

p p oq 

where a = a{t) is a new arbitrary function of time. Equation ([39|) express Co in 
terms of C 2 and V. 


The last equations, corresponding to the linear term in momentum and the 
remaining term, involves Ci and V. After some algebra, we can transform this 
system into an equivalent one, consisting of an equation for C 2 only. 


dC2 ^ „ dC2 

dt ^ dq 


a+-{q-a), 

P 


(40) 


and an equation determining the potential. 


dV dV 
dt ^ dq 


P 


V + {ap + (ap + pa)p/p- 2pa - pa)- 


P 


(- + ^)|-(-(?-^) + ^)C2 

p p^ 2 p 

1 d 


p^ dt 


(pVo(t)), 


(41) 


where Vo(^) is an arbitrary function of time. 


The strategy for solving is clear. First, we have to solve (^D]), obtaining 

C 2 in terms of {q,t). Then, inserting this solution into (^p, we arrive at a well 
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defined partial differential equation for the potential. If its solution is available, 
we can obtain Cq and the corresponding associated compatible vector held using 

(S). 


Fortunately, the solution for equation 
a, and reads 


C2 = - (p(g -a)+p& + Q) , 
P 

where Q = Q{q,t) is implicitly dehned according to 


is available for arbitrary p and 

(42) 



(43) 


where F is an arbitrary function of the indicated argument. As F depends on Q, 
the solution has indeed an implicit character. However, the implicit function the¬ 
orem assures that we can always locally solve for Q{q,t) under suitable conditions 
on F. With Q{q, t), we can write C 2 using (^21) and then the equation ( 0 ) for the 
potential. Note that the appearance of implicit relations is not new in the theory 
of integrable systems (See, e.g., [^, |T7|). 


Taking into account the results of the section, we can write the compatible 
vector held as 

/P _ pidV/dq + a)+ p{q-a) \ ^ 
dq Vp p{p-a)-p{q-a)-Q J dp' 

To obtain explicitly the compatible vector held, we have to dehne the function F in 
(1^ ) and the hnd the potential solving (^). Remembering that the latter is a linear 
partial diherential equation with two independent variables, we conclude that an 
additional arbitrary function will appear after solving (^) by characteristics. In 
the continuation, we illustrate the whole procedure with specihc examples. 


4.1 Quadratic invariants 

Let us hrst consider the case of invariants quadratic in momentum. For this class 
of solutions, we take F = 0 in (^5[). According to (|12D, we then have 

C 2 = & + -{q - a). (45) 

P 

Inserting this form into ( 0 ), we obtain the equation determining the potential. 
The solution is 

V = Vo{t) + {ap-paf--Pf + -^u[^-^\ , (46) 

P 2p P^ \ P ) 
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where U is an arbitrary function of the indicated argument. This is an example 
of how the solution of the equation (^I]) by characteristics yield an additional 
arbitrary function in the complete solution. According to (H), the compatible 
vector held is 


d 


V = 


P 


dU/dq 


d 


dq \p J dp' 


where 


q = 


9 - ^ 
P 


(47) 


(48) 


Equation (|15|) indeed dehnes the class of potentials admitting an invariant 


quadratic in momentum [[1^. To proceed with Frobenius method for the derivation 
of invariants, we use a new momentum variable 


p = p(p-a)- p{q - a ), 
so that the compatible vector held reads 


- ^ (-A _ 

pp \ dq dq dp 


The integral surfaces of v are 


(49) 


(50) 


P 


f = Y + U{q), g = t. 
This yields directly the energy-like invariant 


(51) 


7= Y + ci®. 


(52) 


since, incidentally, df/dg = 0 along the trajectories. But I is the invariant 


quadratic in momentum derived in |]^, expressed in terms of transformed co¬ 
ordinate and momentum. Actually, (|48| - ^9D is an example of generalized canonical 
transformation [^. In the present context, I as given in (|^) corresponds to an 


integral surface of the compatible vector held. 


4.2 Giacomini potentials 

As another particular choice of arbitrary functions, consider 

p=l, a = 0, Vo = 0. (53) 
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Examining equations (p9D and (|iT|), we get 


^0 = - 


dV 


C2 = 


dV/dt 


(54) 


dq ' ^ dVjdq 

Also, instead of using expression for C 2 , here it is more useful to rewrite (|i0| ) 
as 

dC2/dt dV/dt 

dC 2 /dq dV/dq ^ ’ 

The latter result simply states that C 2 and the potential are not functionally 
independent, C 2 = (^ 2 ( 1 ^)• Thus, the potential satisfy 

dV dV 

^4C-,(r)^ = 0. (56) 

According to (|4^) , the compatible vector field can be expressed as 

^ ^ (_d _ 1 

dq p — C 2 {V)dp) 

By inspection, one of the characteristics of v is a function of p and V only, 

f = f{p,V), (58) 

where 


df 


df 


= 0 


dV p — C2{V) dp 
Using the second characteristic g = t, we have, by virtue of (^) and (|58| - ^^ , 

df 


dg 


= 0 , 


(59) 


(60) 


so that f{p, V) is a hrst integral. 


For general C 2 (U), the solution of (|59|) is not known, so that Frobenius method 
only proves weak integrability. However, there is a connection with the work of 
Giacomini [|^, in which it was made a search of invariants depending on the 
momentum and the potential. Giacomini has found exactly equations (^6]) for 
the potential and (|59| ) for the invariant. Giacomini |]^ as well as Bouquet and 
Lewis 10 were able to construct explicit solutions for this system of equations by 
choosing suitable functions G 2 (U). Nevertheless, we have proven the interesting 
result that the Giacomini potentials are always at least weakly integrable, for 
arbitrary G 2 (U), in the sense that they admit a compatible vector field. Also note 
that Giacomini’s solution specified by p = 1, cr = 0 and Vq = 0 does not excludes 
the case of invariants quadratic in momentum, specified by F = 0. However, the 
intersection of the two cases (p = 1, a = 0, Vq = 0 and F = 0) only yields the 
trivial result that energy is conserved for a time-independent potential. We do not 
pursue any further the use of Frobenius method for Giacomini potentials, since it 
will not give new results. 
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4.3 An infinite class of weakly integrable potentials 


As a final example, we consider the choice 

F{s) = sjk 


( 61 ) 


for the arbitrary fnnction in eqnation (|43|). Here, s denotes an arbitrary argnment 
of F and fc is a nonzero nnmerical constant. The choice m allows obtaining Q 
globally nsing (^31), and then C 2 by (|^). The resnlt is 


a = (T 


p \ pj p^ , 


(62) 


It is possible to proceed to the solntion of (|iT|) nsing this fnnction C 2 , for arbitrary 
p, a and Vq. However, the result is awkward and we content ourselves in taking 


a = 0 , Vo = 0 


(63) 


in the continuation. Thus, using (|4lD and (| 6 ^, we hnd the following equation for 
the potential. 


dV q 


A dV 2p 


P , ^PP , 2p 


FF + -AP+-AT + k)-^ ^ = -^V- ^ + ^ + ^{T + k)-^r-, (64) 


dt p 
where 


P 


dq 


P 


T = 


P P 




P dt' 


The solution for 


IS 


r = r(()r + -u(q ), 

p 


(65) 

( 66 ) 


where 


T{t) = 


2p2 


exp 2 


c*' dt"/p^\ ( p 3pp 2p 


T + k I \ p ' p 2 p3 


H- 5 —I—t(T + k) dt' (67) 


and U is an arbitrary function of 


^ q ( [^ dt' /p^' 

q = - exp — 


P 


The associated compatible vector held is 


T + k 


( 68 ) 


V = 


d (p 


pdV/dq + pq 


d 


dq \p pp — pq — {T + k)~^q/p I dp ^ 


(69) 
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with V given in (|6^). 


To try to find the integral snrfaces of v, let 


p = pp - pq - {T + k) ^q/p. 


( 70 ) 


In terms of p and q given in (|68|) , we express the characteristic equation associated 
to V as 


p^ + {T + k) ^ 

dq 


exp 


dt'/p^ \ 
T + kj 


p + p'^p exp 2 


dt'/p^ \ 

T + kj 


q + 2p^Tq + 


dU 

dq 


= 0 . 

(71) 


Remembering that t is simply a parameter in this equation (t = cte. is one of the 
integral surfaces of the compatible vector field), we can identify (^) as an Abel’s 
equation of second type. Such Abel equation of second type is not solvable in 
terms of elementary functions for arbitrary U{q). Hence, in general the class of 
potentials (|66D is only weakly integrable. Finally, we observe that, if p = 1, one 
obtains the Giacomini class of solutions. However, for p ^ 1 the potential (|66|) is 
new. 


5 Conclusion 


In this work, we have presented Frobenius method as an attractive strategy for 
the integrability analysis of one-dimensional time-dependent Hamiltonian systems, 
which may be cast in the form of three-dimensional dynamical systems. We have 
shown that the compatible vector fields to be determined can be considered to 
depend only on a function C{q,p,t). The compatibility condition then leads to 
the basic equation ([T^), determining both C and the potential. Starting from 
some hypothesis on the functional dependence of C, we have been able to derive 
some already known integrable potentials, namely the potential with an invariant 
linear in momentum, Sarlet’s potential, the potential for systems with a quadratic 
invariant and Giacomini’s class of solutions. However, the potential for the system 
with a quadratic invariant and Giacomini potentials are only particular examples 
among a richer class, defined by the solutions of (^T]). These new solutions depends 
on the arbitrary functions p, a, F and a further arbitrary function arising after 
resolution of (^), besides the trivial function Vo(^)- In the final part of Section 
IV, we have constructed a particular weakly integrable system contained in this 
family. 

Frobenius method only assures weak integrability, without providing a first in¬ 
tegral in all cases. In this respect, there is a similarity with Painleve analysis Q, 
which only points out when integrability is likely to occur. Indeed, in Painleve 
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analysis the construction of an invariant is a separate issue, to be addressed after 
the singularity structure of the solutions is understood. Also note that, once a 
compatible vector held is found, the Frobenius procedure for construction of in¬ 
variants is entirely analogous to Lie’s ii approach. The difference between the 
two methods is that a compatible vector held is not necessarily a generator of Lie 
symmetries. 


Even if the existence of a compatible vector held does not allows the construction 
of a hrst integral, if a hrst integral is known it is given by Frobenius method. 
Indeed, let us make the transformation 


dJ/dq 
dJ/dp 


(72) 


on equation ([I^), where J = J{q,p,t) is an arbitrary function of p, q and t such 
that dJ/dp 7 ^ 0. A little calculation using the form (17^) shows that 

“(C) + ■ (73) 


Hence, if J is an invariant, then the right hand side of (^) is identically zero and 
the basic equation (pTf) is satished by the ansatz (^). Moreover, the choice ([7^) 
readily implies that v( J) = 0, showing that the hrst integral J is also a hrst integral 
of the compatible vector held. This shows that all invariants for one-dimensional 
motion under a time-dependent potential may be furnished by Frobenius method. 
In this sense, Frobenius method has an universal character. 


In most situations, no hrst integral is available a priori. In these cases, the 
crucial point for the ehectiveness of the Frobenius approach is a judicious choice of 
the form of the compatible vector held. Clearly, it is possible that other functional 
dependence of the coefficient C on the compatible vector held, diherent from the 
choices made in the present work, may lead to useful results. Another question 
that deserves attention is the extension of Frobenius method and the notion of 
weak integrability to higher dimensions. To conclude, we have seen that Frobenius 
method is a powerful tool in integrability analysis, and we expect that further 
results can be produced by its systematic application. 
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